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Abstract. Wc find discrete holomorphic parafermions of the Ashkin-Teller 
model on the critical line, by mapping appropriate interfaces of the model to 
the 0(n = 1) model. We give support to the conjecture that the curve created 
by the insertion of parafermionic operators at two points on the boundary is 
SLE(4, p, p), where p varies along the critical line. 



PACS numbers: 05.50.+q, 11.25.Hf 

1. Introduction 

The discovery of the Schramm-Loewner evolution (SLE) by Schramm [1] opened a 
rigorous way to study conformally invariant systems. Some weh-known results from 
physics literature were proved [2, 3, 4], new formulas were also discovered [5]. Notably, 
conformal invariance was proved for percolation clusters [2] and for Ising Fortuin- 
Kasteleyn and spin clusters [3, 4]. One of the key steps in these proofs is to find an 
appropriate discrete holomorphic parafermion for the critical interface. Relating the 
holomorphic parafermion to a specific property of the critical curve, one can conjecture 
and then prove that the continuum limit of the interface is SLE. Different methods 
of finding discrete holomorphic parafermions for statistical models and their relation 
to integrability were investigated by Cardy and collaborators in a series of papers 
[6, 7, 8]. In some cases [6, 8] the relation of the parafermion to a critical interface is 
known, and so there are conjectures on the SLE corresponding to the continuum limit 
of the interfaces, but in other cases [7] the corresponding interfaces are not known and 
so there is no known SLE. 

One of the interesting models where neither discrete holomorphic parafermions 
nor corresponding SLE were known up to now is the Ashkin-Teller (AT) model [9, 10] 
(the phase diagram is described in [11, 12, 13]). From the Conformal Field Theory 
(CFT) point of view, the AT model is interesting because it has a critical line with 
constant central charge c = 1 and changing critical exponents. Although it is possible 
to define different kinds of critical interfaces for the AT model [14, 15], it is shown 
numerically [14] that some of the natural possibilities are not related to simple SLE's. 
Recently, in [16], the AT model was studied on iso-radial graphs, the critical surfaces 
on generic iso-radial graphs were found, and discrete holomorphic parafermions were 
defined at some particular points of the critical line, by using algebraic relations 
between spin and disorder variables. In the present paper, we introduce discrete 
holomorphic parafermions all over the critical line, without using these relations 
directly. The idea is to map some particular interfaces to the 0(n = 1) model on 
the square lattice, and then to exploit some of the results of [8] . 

Our main conjecture is based on the continuum Gaussian theory for the Solid-On- 
Solid (SOS) model associated to the AT model, and on the known relation between 
this theory and SLE4 [17]: in the upper half-plane, we expect the curve created by 
the insertion of a parafermion at the origin and at infinity to have the statistics of 
SLE(4, — 1, y/g — 1), where g is the coupling constant of the Gaussian theory (see 
Section 4 for details on the definition of g). 

The structure of the paper is as follows. In Section 2, we recall the definition 
of the AT model on the square lattice and its mapping to the staggered six-vertex 
model. In Section 3, the six-vertex model is mapped to the 0(n ~ 1) model, and, 
using this mapping, we find a lattice holomorphic parafermion for the AT model all 
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over the critical integrablc surface. In Section 4, we formulate conjectures on the 
relation of our interfaces to SLE. Section 5 is dedicated to bulk critical exponents 
in the 0(n = 1) loop model associated to the AT model. The results are checked 
numerically by transfer-matrix diagonalisation. 

2. The Ashkin-Teller model 

2.1. Definition and graphical expansion 

The Ashkin-Teller model can be defined on any graph, but, for clarity, we will restrict 
to the square lattice in this paper. At each vertex j of the lattice, we define two spin 
variables Cj and Tj, which can take individually the values ±1. A spin configuration 
{(Tj,Tj} gets the Boltzmann weight: 

Y[ W{i,j) , where W{iJ) = exp [j3^j{a^aj + TiTj) + a.,j a^ajnTj] , (I) 

and aij,l3ij = ax,Px (resp. ay,l3y) if is a horizontal (resp. vertical) bond. 

We denote by (...) the averaging with respect to the normalized Boltzmann weights 
W{i,j)/Z, where 

^= E n^(*'j)- (2) 

Let us recall the graphical expansion of this partition function [18, 19]. Introducing 
the change of variables {<Jj,Tj) — >■ ((7j,Tj), where rj = <JjTj, the edge interaction can 
be written as: 

W{i,j) = e"-^.'^.' [cosh/3,,(l + T^rj) + a,a, sinh + r/rj)] . (3) 

First, we fix the values of the rj spins, which define a domain-wall (DW) configuration 
on the edges of the dual lattice. Let {i,j) be an edge of the original lattice. If (i, 
is crossed by a DW, then it gets a weight e~"'J . If not, then it gets the weight 
e"'J [cosh2/3ij -t- aicyj sinh2/3y ]. We depict the first term as an empty edge, and the 
second one as an occupied edge. Thus the partition function reads 

Z = const X ^ (tanh2/3:r)^"^^'(tanh2/3j,)^''('^' 

G\G' 

X (e^"" cosh2/3y)-^==('=^')(e2"- cosh 2/3^ )-^''(^') , (4) 

where G (resp. G') is a subgraph of the original (resp. dual) lattice, all vertices 
in G and G" must have even degree, G and G' must not intersect each other, and 
£x{G),£y{G) are the numbers of horizontal and vertical bonds of G. An example 
configuration is shown in Figure 1. 

2.2. Mapping to the six-vertex model 

One interest of this graphical expansion is that it maps to the six- vertex model on the 
medial lattice! . The correspondence is given in Figure 2. The Ashkin-Teller model 

X An earlier construction, described in [12], relates the AT model to a staggered eight- vertex model, 
without using the r' variables. However, in the present paper, we focus on the mapping to the 
six-vertex model because the latter has a known discrete holomorphic parafermion. 
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Figure 1. An example configuration in the graph expansion (4). The spin 
variables live on the vertices of the underlying square lattice. Domain walls for 
the rj spins are depicted in thick, black lines, and high-temperature polygons 
for the cTj variables, in thick, grey lines. Crossings between the black and grey 
clusters are forbidden. 



is critical when the six-vertex model is not staggered. The weights of the six-vertex 
model with parameter A — (a^ + ~ c^)/{2ab) — — cos A can be paramcterised as 

a, fe, c = sin(A — u), sinu, sin A . (5) 

This corresponds to the weights of the AT model: 

, „^ sinu e~^""= sin(A — w) 

tanh 2/3^ ^— , — — — \ ^ ' (6) 

sin A cosh2p2; sin A 

, „^ sin(A — u) e~^"f sinw 

tanh2/3.„ = — \ ^ ' , — — = ^— . 7) 

smA cosh2py smA 

The variable m is a spectral parameter, whereas A determines the universality class. 
The isotropic point is at u = A/2. Some special values are: Ap = 0, Apz = f j -^i = f 
and A = corresponding to the four-state Potts, Z4 Fateev-Zamolodchikov [20], 
Isingxising and XY models. When A is varied, the central charge remains constant 
c = 1 , but the critical exponents change. For example, the correlation exponent is 

, , 27r - 2A 

More generally, it was shown in Ref [16] that the critical weights of the AT models 
on any Baxter lattice are given by (6)-(7), with the spectral parameter u related to 
the angle of 9 the rhombic faces by 9 — ttu/X. In the next Section, we will describe a 
discrete holomorphic parafermion in the critical square-lattice AT model. Note that 
this parafermion is also present for the AT model on a Baxter lattice, at the critical 
value of the Boltzmann weights. 

3. Discrete holomorphic parafermion 

In [16], using the spin variables of the AT model, discrete holomorphic parafermions 
were found at some particular points of the critical line, namely at the four-state Potts, 
FZ, Ising and XY points. In this Section, we will exhibit anot/ier discrete holomorphic 
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Figure 2. Mapping from the graph expansion of the AT model to the six- vertex 
modeh The colour of the possible graph bond depends on the sublattice of the 
medial lattice. 





Figure 3. Arrow configuration corresponding to the graph configuration shown 
in Figure 1. The arrows of six- vertex model live on the edges of the medial lattice. 



parafermion all along the critical line of the AT model. This parafermion is defined in 
the loop formulation of the AT model. The key idea is to use the chain of mappings: 

Ashkin- Teller — > Six- vertex — > 0(n = 1) . 

The first mapping is described in Section 2, and the second one is given, for example, 
in [8]. For completeness, we recall it in Figure 4. The resulting loop model has seven 
possible vertices, and each closed loop has a Boltzmann weight n = 1. We will denote 
this model the 0{n = 1) model for short. 

It was shown in [8] that the 0(n = 1) model, with weights given by (5), possesses 
a discrete holomorphic parafermion ip{z)^ associated to the insertion of one path at 
point z. The spin s of this parafermion is related to the parameter A by: 



In terms of the AT model, is defined on the edges of the covering lattice (the 
union of the original and dual lattices). If z denotes an edge of this lattice, the operator 
tjj{z) inserts a a variable and a r' DW at the ends of the edge z. The two-point function 
{il^{zi)^(z2)) also contains a non-local phase factor given by the simultaneous winding 
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Figure 4. Mapping between the six-vertex model and an 0(n = 1) loop model. 
The correspondence between the arrow and loop vertices shown here is valid on the 
sublattice 1 of the medial lattice. On the sublattice 2, the same correspondence 
holds, with all arrows reversed. 



of the r' DW and cr cluster between zi and Z2. We might decompose symbolically 
into three factors: 

■0 ^ e*^ xax ^r' , (10) 

where is a contribution to the winding angle (see above) and fir' creates a r' domain 
wall. 

Let us interpret the effect of -ip on the AT model, when inserted at two boundary 
points a and b. In this situation, there must be. in the 0(n = 1) model, a path 7 
going from a to 6 (see Figure 5). The occupied area adjacent to 7 is bounded on 
one side by a r' domain wall, and on the other side by a cr high-temperature cluster. 
So {ip{a)ip{b)) corresponds to the two-point function (a{a)a{b))i,ab, where (. ..)i,a6 
denotes the averaging with Boltzmann weights W{i,j)/Z and the following boundary 
conditions (denoted BCi): free boundary conditions for cr and t' = 1 on one arc (ab) 
of the boundary, r' = — 1 on the other arc. 

We can also think of boundary conditions which allow the path 7 to end anywhere 
in a given interval of the boundary. Let a, 6, c be three marked points on the boundary. 
We consider the two-point function (cr(a)cr(6))2,a6c, where (. . ■)2,abc is the Boltzmann 
average with the following boundary conditions: free BC for cr, r' ~ +1 on [ab) and 
(6c), r' = —1 on (ac). In the 0(n = 1) model, this forces a path 7 to go from a 
to a point on (6c), as shown in Figure 5. We call these boundary conditions BC2. 
If /i^' is viewed as a boundary-condition changing operator, the correlation function 
(cr(a)cr(6))2 can be written as a three-point function, reflecting explicitly the role of 
point c: 

{a{a)a{b))2,abc = (^(a)a(6)M,. (c)) (11) 
4. Relation to SLE 

The six-vertex model can be transformed into an SOS model by introducing height 
variables ip on the faces of the medial lattice, so that two neighbouring (/j's differ by 
±7r/2, the highest being on the left of each arrow. The SOS model rcnormaliscs to a 
Gaussian theory with action 

.4 = A 1 1^,1 v.. » = ^^. (12) 

When the AT model is defined on a system with one or more periodic directions 
{e.g., a torus or a cylinder), the height ip is only well-defined up to the identifications 
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Figure 5. Random curves associated to the parafcrmion ili, defined on a domain 
f! with boundary conditions BCi (left) and BC2 (right). In both cases, closed 
loops may be present around 7. 



if = if + TT,ip = —(p. The first corresponds to a height defect around one direction of 
the system, while the second is a 'twist' induced by an odd number of DW's winding 
around the system [19]. The theory with action (12) and the above identifications 
is called the Z2-orbifold of the compact boson. However, in this Section, we will 
be dealing only with simply-connected domains, so the SOS configurations are always 
well-defined without any identification, and the continuum limit is simply the Gaussian 
model (12). 

In the case of boundary conditions BCi, the presence of the path 7 induces a 
height gap at a and b. So we expect the SOS model to be a free field with Dirichlet 
BC, and boundary values — fi on one arc, ip2 on the other arc. Since the DW and high- 
T cluster play the same role at the critical point, we must have fi ~ ip2- Schramm 
and Sheffield have shown [17] that the contour line in such a model is SLE(4, pi, P2), 
where pi = ^p\l^>* — l,/92 = '^2/'^* ~ 1, and (/?* is a universal constant. Now, we 
use the results of [27, 21]: when a consistent normalisation is chosen for we have 
tp* = 11 / \/Ag and the operator inserting the height gap 5lp has conformal dimension 




Comparing to the spin (9) of the parafermion s = .g/4, we obtain ip^i = (p2 ^ 7r/2, and 
so we conjecture that the curve 7 has the statistics of SLE(4, ^ — 1, ^ — 1). 



5. Critical exponents and fractal dimension 

In this Section, we derive the bulk exponents for watermelon correlation functions in 
the 0(n = 1) model, and give numerical results on some of these exponents. The 
0{n = 1) model is a peculiar loop model, where the central charge is fixed, but the 
critical exponents vary along the critical line. 

To calculate bulk exponents, it is easiest to consider the loop model on a cylinder 
of circumference L sites {L even). In this setting, the model is described by the Z2- 
orbifold theory, and part of the conformal spectrum [19] is given by 'electromagnetic' 
exponents 

2 2 
e am 

zg 2 
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In the transfer-matrix formalism, the fc-leg watermelon exponent corresponds to 
the dominant eigenvalue in the sector with k strands propagating along the cylinder. 
For a given k, there are several exponents, according to the parity of the sites where the 
k strands sit. Indeed, because of the staggering in the 6V/0(n 1) correspondence, 
strands sitting on an even (resp. odd) edge are oriented positively (resp. negatively). 
Each strand contributes a height defect of magnetic charge m = ±1/2, so if we write 
fc = fci +/c2; where fci, k2 are the numbers of strands sitting on even and odd edges, the 
total magnetic charge is m = (fci — fc2)/2. Hence, when fci ^ k2, we get the exponent: 

g(fci - fca)^ L / f, nA\ 

Xki.k2 = ^ , ki^k2. (14) 

In the special case fci — k2, the total magnetic charge vanishes. Our numerical results 
allow us to conjecture to exponent: 

Xij = Y . (15) 

The fractal dimension of the path 7 is related to a two-leg watermelon exponent: 
df = 2 — X. For fc = 2, there are two choices: X2.0 = g/2 and i ~ 1/2. The choice 
Xi^i gives a fractal dimension df = 3/2, which is the correct value for SLE4. Note 
that the direct numerical calculations of Picco and Santachiara [23] on the AT model 
show that the fractal dimension of the boundaries of r' clusters is also 3/2. 

The central charge and the exponents of the loop model are extracted using the 
finite-size formulae for the eigenvalues Aj of the transfer matrix [22]: 

- log A, ^ i/^ + ^ (-^ + X,) . (16) 

In Figures 6-12, we depicted numerical data for the central charge, the thermal 
exponent and some watermelon exponents of the 0{n — 1) loop model. 



6. Conclusion 



In this paper we first found a discrete holomorphic parafermion for AT model on 
the whole critical line. There may exist other discrete holomorphic parafermions in 
this model, but the one we describe is defined in terms of the 0(71 = 1) loop model, 
which enables us to relate it to the SLE model. Our conjecture is compatible with our 
numerical calculations on the 0{n = 1) model. Of course, more precise numerical 
calculation is needed to confirm our conjecture, such as the left-right Schramm's 
formula. Unfortunately, no Schramm's formula for SLE(k, pi, ^2) is known, but for 
K = 4 it could be tractable, thanks to the relation with the Gaussian Free Field. One 
can also think of Monte-Carlo simulation as it was done for the Ising model [25], and 
martingale arguments [26]. 

At the point A = 7r/2, the AT model is equivalent to two decoupled Ising models. 
Oshikawa and Affleck have studied the boundary CFT of this model, in relation to the 
Z2 orbifold theory [27]. In particular, our parafermionic observable -0, when inserted 
on the boundary, corresponds to a jump in Dirichlet boundary conditions for the 
orbifold {(fa = — tt in the notations of [27]), together with the insertion of a cr 
operator. An interesting direction would be to try and extend these results to the 
critical line of the AT model, and analyse the null-vector equations of the boundary 
operators, to relate them properly to SLE. 
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Figure 6. Central charge of the loop model. The expected value is c = 1 on 
whole critical line. 




Figure 7. First and third thermal exponents of the loop model. There is a 
crossing between the exponents X = 2 ~ l/i'{X) and X = 1. 
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Figure 12. Six-leg watermelon exponent X^^s of the loop model. 
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